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The anomalous forward-backward asymmetry (Afb) in tt production at the Tevatron is in ap- 
parent contradiction with the tt charge asymmetry (Ac) measurements at the LHC, which agree 
well with the standard model predictions. We argue that associated production of a state with [ut\ 
flavor quantum numbers can lead to a sizable negative contribution to Ac- Exchange of such a 
resonance in the t— channel would lead to positive contributions to Afb and Ac, as has been exten- 
sively discussed in the literature. Given the additional negative Ac contribution, both the Tevatron 
and LHC data can be naturally accommodated within this framework. A simple realization of this 
setup is the well known example of a Z' with flavor off-diagonal up-top couplings. We provide a 
detailed study of this model, demonstrating that it indeed reproduces the tt asymmetry data, and 
is compatible with other constraints, e.g. t ± jet resonance searches, tt inclusive jet multiplicities, 
and atomic parity violation. 



I. INTRODUCTION 

The inclusive forward-backward asymmetry (Afb) 
in the tt rest frame has been measured by both the 
CDF [1, 2] and D0 [3] collaborations and found to be 
significantly larger than the standard model (SM) predic- 
tion. The naive average in semileptonic tt events, adding 
the uncertainties in quadrature, is 

A FB = 0.174 ±0.038. (1) 

The NLO QCD predictions, including leading elec- 
troweak (EW) contributions [4-7], are [7] 

Af£ = 0.088 ±0.006, (2) 

when using the LO total cross section in the denominator 
of the asymmetry. The error is due to scale variation with 
the central value evaluated at \i = mt and the low and 
high values evaluated at /i = mt /2 and 2m t , respectively. 
Both CDF and D0 have also measured the asymmetry 
in bins of the tt invariant mass and tt rapidity differ- 
ence. Only CDF [1, 2], however, unfolds to the partonic 
("truth") level obtaining 

A l p% = A FB (m t t < 450 GeV) = 0.078 ± 0.054 , (3a) 

A^f = A FB (m t t > 450 GeV) = 0.296 ± 0.067 , (3b) 

to be compared with the SM (NLO QCD and EW) pre- 
dictions [4, 5, 7] 

(A^) SM = 0.0621°;™ (4t) SM = 0.129t° Z, (4) 
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again using the LO cross sections in the asymmetry de- 
nominators, with the errors due to scale variation as 
above. At the top reconstruction level, D0 does not 
observe a significant rise in the asymmetry vs. the tt 
invariant mass, however an increase vs. the tt rapidity 
difference is observed [3]. 

A related observable at the LHC is the charge asym- 
metry in tt production, Ac- In contrast to the forward- 
backward asymmetries, the measurements of Ac agree 
with the SM expectations. For example, the average of 
the inclusive ATLAS [8] and CMS [9] measurements at 
= 7 TeV, 

A c = (1.15 ±1.25) • HT 2 , (5) 

agrees within errors with the SM prediction [5, 7] 

A% M = (1.23 ±0.05) • 10~ 2 . (6) 

ATLAS collaboration also presented the first measure- 
ments of Ac in bins of m t i [8] 

A l g w ee A c (m t t < 450 GeV) = -0.053 ± 0.088 , (7a) 
A^ gh = A c (m t i > 450 GeV) = -0.008 ± 0.047 , (7b) 

which are consistent with the SM predictions [5] , A)% w = 
0.015 ± 0.001 and A^ gh = 0.026 ± 0.001 (all LHC asym- 
metry predictions correspond to LO cross sections in the 
denominators, and the errors are due to scale variation, 
as above). Compatible results, but in three m t i bins, 
were presented by the CMS collaboration [9]. 

The discrepancy between the SM predictions and the 
measured asymmetries at the Tevatron could be due to 
an unknown QCD effect, or an unidentified experimen- 
tal bias. Alternatively, it might be a hint of dynamics 
beyond the SM (for a review, see, e.g., [10]). However, 
new physics models are challenged by the requirement of 
achieving simultaneous agreement with the anomalously 



large Afb measurements and the SM-like Ac measure- 
ments, as the underlying candidate NP processes often 
yield correlated positive contributions to both [12-14]. 1 

The main purpose of this paper is to point out that 
there is a class of models in which positive contributions 
to Afb and Ac are correlated with an additional nega- 
tive contribution to Aq. Therefore, the latter tends to 
be small, in agreement with current data. A state with 
[iu] or [id] flavor quantum numbers is required. The pos- 
sible examples are the [iu] flavored Z' [17] (see also [18]) 
or [id] W [19] vector mediators, and the [iu] flavored 
SU(2)l doublet scalar mediator (j) [20]. In this paper 
we focus on the Z' color- and weak-singlet, with a cou- 
pling to the right-handed up and top quarks. Exchange 
of the Z' in the ^-channel, see Fig. la, would lead to an 
increase in vs. m t t, and a positive contribution to 
Ac in excess of the measurements [12, 14], due to for- 
ward peaking. However, associated production of the Z' 
with a top-quark would produce an additional negative 
contribution to Ac- Presumably, this would also be the 
case for the W and <ft. However, we leave the question 
of the overall viability of these models for future studies. 

The leading order ug — >• Z't — » uit Feynman diagrams 
are shown in Figs, lb and lc. The effect we are interested 
in is due to the top quark exchange diagram in Fig. lb. 
The Z' decay yields a i quark which tends to be boosted 
in the same direction as the incoming u quark. Taking 
into account the harder u quark vs. gluon PDF's in the 
proton, one concludes that on average the i is produced 
with a larger rapidity than the t, thus yielding a negative 
contribution to Ac- The efficiency of this mechanism is 
illustrated in Fig. 2 for one of the Z' benchmark points 
(p 3 ) listed in Table II. The pp (Z't, Z'^i) itX differ- 
ential cross section plots in Fig. 2a and Fig. 2b, respec- 
tively, exhibit the dominance of i (t) production at larger 
(smaller) rapidities and the dominance of ti production 
for \y t \ — \yt\ < 0. According to Fig. 2b, the charge 
asymmetry from associated vector mediator production 
alone can be large and negative, e.g., Ac — —18% for the 
example shown. 

At the LHC, the cross section for the CP conjugate 
process, ug — > Z'H — >> uti, is typically an order of mag- 
nitude smaller, due to the iZ-quark PDF in the initial 
state. This asymmetry has been studied previously in 
[21]. A corollary is that the negative contribution to Ac 
does not depend on whether the Z' is self conjugate or 
not. However, bounds on same sign top production can 
rule out a self conjugate Z' . In fact, flavor symmetric 
realizations of the vector ^-channel models [22, 24, 26] 
proposed to trivially evade the bounds on same sign top 
or single top production and FCNC's, e.g., D — D mix- 



However, in the case of a light s-channel axigluon, it has recently 
been shown that different couplings to the u and d quarks can 
lead to partial cancelations between the uu and dd contributions, 
thus sufficiently suppressing Ac, with only a small impact on 
Afb, which can be sizable and positive [15, 16]. 
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Figure 1: Feynman diagrams for (a) t-channel Z' exchange 
contribution to uu — »• ££, and associated single Z' production 
in the £ channel (b) and s channel (c). 



ing, would contain a CP conjugate pair of [iu] and [tu] 
flavored Z n s. At the Tevatron, associated production of 
the vector mediators produces a negative contribution to 
Afb- However, this effect is suppressed relative to the 
positive Afb contribution from Fig. la by the smaller 
gluon vs. ^-quark PDF's inside the proton. 

It is worth pointing out that in models for enhanced 
Afb in which a mediator is exchanged in the ^-channel 
rather than the t-channel, e.g., a color triplet or sextet 
scalar S [24, 26-29] with couplings StRUR (where the 
color indices have been suppressed), single mediator pro- 
duction via the analog of Fig. lb would lead to a positive 
contribution to Ac [30]. This is easily understood, as in 
this case the mediator would decay to a top quark rather 
than an anti-top quark. 

Associated Z' production is subject to several con- 
straints which need to be checked in order to estab- 
lish the viability of our mechanism for reducing Ac- 
The most relevant ones are: i) the LHC ti cross sec- 
tion measurement, which constrains the product <j(pp —> 
Z't) x Br(Z' tu) [31]; ii) the CDF, CMS and ATLAS 
collaborations Z' — >• t + jet resonance searches [32-34], 
that yield direct bounds on the above product [35]; Hi) 
the CMS measurement of the jet multiplicity distribu- 
tion in semileptonic ti events, which is consistent with 
SM Monte Carlo studies [36] , thus potentially constrain- 
ing NP models which produce extra jets [37]. 

Constraints ii) and Hi) become weaker as the Z' mass 
approaches m t from above. Because the u quark be- 
comes softer, the resonance searches lose sensitivity and 
the jet multiplicities become more SM like. Furthermore, 
a reduction of Br(Z' — >• tu) trivially reduces the impact 
of searches i)-iii). In fact, we will see that the ranges 
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Figure 2: The relative differential cross sections (a) 
1/a da/d\y ti t\ and (b) 1/a da/d(\y t \ - \yt\) for pp -> Z't(t) 
ttX at ^/s — 7 TeV, for benchmark point p3 in Table II. The 
subscripts B and F refer to \y t \ — \yi\ < and > 0, respec- 
tively. 



M z > ~ 200 - 300 GeV and Br(Z' -> tu) ~ 1/4 - 1/3 
are preferred, with the best overall fits occurring for Mz> 
near 220 GeV. A sufficient negative shift in Ac can be 
obtained, after taking the branching ratio into account, 
due to the sizable negative charge asymmetry intrinsic to 
associated Z >• iu) production. 

The favored range for Br(Z' — >• iu) raises an immedi- 
ate question: what are the viable candidates for the miss- 
ing dominant Z' decay? Some possibilities that come to 
mind are invisible decays, decays to light quark pairs, 
and decays to a pair of scalars which subsequently decay 
to quark pairs, which could include a virtual top quark. 2 
The first option yields mono-top events [39, 40]. The 
CDF collaboration performed a search for the monotop 
signal and obtains a bound on the cross section of 0.5 pb 
based on 7.7 fb -1 of data [41] (for a Z' mass between 
and 150 GeV, while the bound would be weaker for a 
heavier Z'). This means that all of the benchmarks to be 
discussed below are allowed to decay predominantly into 
invisible states. The second and third options result in a 
t + n jets final state (n > 2). Single top searches at the 



2 The last option could be realized in flavor symmetric models 
based on new strong interactions, in which the Z n s are composite 
elements of a flavor nonet, analogous to the K*± of QCD, with 
dominant decays to pseudo Nambu-Goldstone bosons [38]. 



Tevatron [42-45] and the LHC [46, 47] have not been op- 
timized for the higher number of high pt jets compared 
to the SM. Nevertheless they may pose an important con- 
straint. A detailed investigation of these issues is left for 
future work. 

Finally, the Z' model is subject to low energy con- 
straints from atomic parity violation (APV) measure- 
ments [54], as are the W and <p models. As it stands, the 
Z' model only contains massive vectors with flavor vio- 
lating couplings to quarks. Thus, a UV completion would 
be required in order to render it renormalizable, and sim- 
ilarly for the W model. As we shall see, at the level of 
an effective theory in which one only considers the effects 
of the Z' ', the severity of these constraints significantly 
weaken (we also comment below on possible UV comple- 
tion of our model) . A shift in the Cs weak nuclear charge 
Qw( 133 Cs) due to the nuclear matrix element corrections 
recently discussed in [55] also ameliorate the APV con- 
straints. Thus, one can not claim any tension with the 
APV data in the preferred region of parameter space for 
the Z' model. 

The organization of the paper is as follows. In Sec- 
tion II we discuss our fits to the Tevatron and LHC ti 
production data, thus arriving at the preferred region of 
parameter space for the Z' model. We will demonstrate 
that the addition of associated Z' production is crucial 
for obtaining a good fit, due to the LHC measurements of 
Ac- This is illustrated in Fig. 7, which is the main result 
of this paper. We also check that the predicted Tevatron 
and LHC m t i spectra agree with the measured ones. In 
Section III we discuss the impact of the tj resonance 
searches and jet multiplicity measurements. Atomic par- 
ity violation is discussed in Section IV. In Section V we 
briefly discuss the implications of the Z' model for the 
leptonic forward-backward and charge asymmetries, and 
conclude in SectionVI. 



II. IMPORTANCE OF ASSOCIATED Z' 
PRODUCTION 

The Lagrangian for the Z' model is given by 

C = g ut u R ^ t R + h.c. + M% Z'J Z'f. (8) 

As previously mentioned, the Z' would not be self- 
conjugate if it transformed non-trivially under flavor 
symmetries introduced to suppress same sign top or sin- 
gle top production and FCNC's [22, 24, 26]. However, 
a self-conjugate Z' would not modify the results of our 
analyses for ti observables. We employ a x 2 function to 
search for optimal ranges of the coupling g ut , the vector 
mass Mz't and Br(Z' — >> tu). Six ti production observ- 
ables enter our fit: the total Tevatron and LHC cross 
sections, <r™ v and cr^ c , the Tevatron inclusive forward- 
backward asymmetry Afb, the ones for m t i > 450 GeV 
and m t f < 450 GeV denoted as A 1 ^ and A 1 ^^ 1 respec- 
tively, and the LHC inclusive charge asymmetry Ac. The 
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Experiment 



SM Theory 



crt T o? V [pb] 7.50 ± 0.48 [60] cr t T f v [pb] 7.067 ± 0.26 [58] 



A FB 

/\high 
^FB 



'%} 17.4 ± 3.8 [2, 3] Aa SM [pb] 0.71,0.48,0.33,0.25 
%] 29.6±6.7[2] Aa SM [pb] 0.46,0.27,0.18,0.13 
%] 7.8 ±5.4 [2] Aa SM [pb] 0.27,0.22,0.15,0.11 



LHC 
^tot 



^t L ot C tpb] 162.6 ±17.1 [59] 



[pb] 172.8 ± 8.2 [61] 
A c [%] 1.15 ± 1.25[8, 9] Aa SM [pb] 2.51,1.60,1.16,0.96 

Table I: Experimental and SM predicted values for tt cross- 
sections and asymmetries at the Tevatron and LHC. Multiple 
entries correspond to Acr SM cross section differences evaluated 
at a = rat/2, rat, 2ra*, 4rat, see text. 



experimental values are listed in the second column of Ta- 
ble I. Where more than one measurement is indicated we 
take the naive average, with errors added in quadrature. 

The NP contributions to the above observables are 
evaluated at LO, whereas the SM contributions are taken 
at the highest order available. The total cross sections 
are expressed as 



Gtot 



tt,z' 



x Br(Z ; 



tu) 



SM 



(9) 



where a tt,z denotes contributions from t-channel Z' ex- 
change, including interference with the SM amplitude at 
LO, and a z 1 denotes the sum of the cross sections for 
Z't and Z'^t production. The inclusive tt asymmetries 
are expressed as 



Act 



tt,z' 



• Aa 



z't 



Aa 



SM 



tt,Z' . Z't w "D ( ryi . T \ i SM,LO ' 



For A FB , 



Aa 1 



.SM 



SM 



SM 

-a B 



A^tt.z' tt,z' tt,z' 
Aa = a F -a B , 

7>\ 



(10) 



(11) 



Aa zt = a F (pp^ Z't/Z'U^ttX) - (a F a B ) ,(12) 

where the F and B subscripts denote events in which the 
top quark is emitted in the forward and backward hemi- 
spheres, respectively, with respect to the proton beam 
direction. Note that Aa z 1 scales as Br(Z' — >• tu). The 
expressions for A 1 ^ and A BB correspond to the restric- 
tions m t i < 450 GeV and > 450 GeV, respectively, for 
the cross sections entering (9)-(12). For Ac, substitute 
pp pp, F \y t \ - \yt\ > and B -> \yt\ - \yt\ < in 
(11), (12), where y t and yi are the t and t rapidities. 
The tt cross section for the z'th m t i bin is written as 



tt,z' 



SM 

(J A 



Gi (pp^ z't/z'U-> ttx) 

(13) 

where the e$ are acceptance corrections, relative to the 
SM ones, which account for a significant decrease in the 
CDF detector's sensitivity to top quark production in the 
forward region. As the forward region features promi- 
nently at large invariant masses in the Z r model, it is 



scale choice 



(m z >,gut) 
[Bi(Z f -^ut) = 



(m Z ',g u t) 
1/4] [Br(Z' -+ut) = 1/3] 



a = rat/2 
a = m t 

ji = 2m t 

a = 4ra t 



(223,0.70) 
pi (221,0.66) 
p 2 (323,0.84) 
p 3 (222,0.63) 
p 4 (372,0.91) 
(224,0.61) 
(374, 0.88) 



(225,0.71) 
(229,0.68) 

p 5 (276,0.74) 

(286,0.73) 



Table II: Benchmark points in (mz f ,g u t) which minimize the 
X 2 for different choices of \i and Br(Z x — > tu), indicated in 
Figs. 3-5 with black dots. If two entries are included, the 
second one corresponds to a higher m Z ' benchmark that is 
marked with a yellow cross in Figs. 4, 5. Benchmark points 
pi, ..,ps (bold) are used in subsequent sections. 



crucial to take these corrections into account when com- 
paring the predicted m t i spectra and the unfolded CDF 
spectrum [22-24]. The t{ are evaluated using the table 
provided in [24] for the approximate CDF tt pair detec- 
tion efficiencies, given in bins of m t t and the rapidity 
difference, Ay = y t — y^. Inclusion of constraints from 
the CDF m t t spectrum would slow down our y 2 search 
procedure substantially, due to evaluation of the e$. For- 
tunately, this turns out to be unnecessary, as the regions 
of parameter space favored by our current search are con- 
sistent with the CDF m t i spectra, see examples below. 
This is also the case for the predicted LHC m t i spectra. 
However, we will see that the LHC is beginning to probe 
the high m t i tail expected in the Z' model. 

In our fits, the NP tt cross sections, and the LO 
SM cross sections entering the denominators of the tt 
asymmetries, are evaluated for a wide range of renor- 
malization scales, \i = m t x (1/2,1,2,4), in order to 
capture the large theoretical uncertainty associated with 
scale variation at LO. The corresponding associated Z' 
production cross sections are evaluated at the scales, 
[i = (mt + m Z ') x (1/4,1/2,1,2). The SM cross sec- 
tions are listed in the third column of Table I. The to- 
tal SM cross section at the Tevatron, cr t s ^ (TEV), is 
identified with the NNLO + NNLL evaluation in [58] 
for m t = 173.3 GeV, and the total SM cross section at 
the LHC, af^ (LHC), is identified with the approximate 
NNLO + NNLL evaluation in [59] for ^fs = 7 TeV and 
m t = 173.3 GeV. At this order, the errors associated 
with scale variation are relatively small. Thus, in our 
fits they are absorbed into the overall theoretical uncer- 
tainty on the SM cross sections. The A<j sm, s entering 
the tt asymmetries are evaluated at the same renormal- 
ization scales as the NP and LO SM cross sections. For 
/jl = mt/2,m t ,2m t they are identified with the values 
given in [7] for m t = 173.1 GeV, corresponding to the 
asymmetries given in Eqs. (2), (4) and (6). They are 
obtained at NLO including electroweak corrections. In 
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Figure 3: % 2 analysis for \i — mt/2, Br(Z ; — > ut) = 1, 1/3, 1/4 
giving la and 2a preferred regions (red), see text. The gray 
area represents the region which is not excluded by ATLAS 
searches for top+jet resonances [34]. The black dot denotes 
the best fit point. 



Figure 4: Same as Fig. 3 for \i — m t , with the yellow cross 
corresponding to the higher mz' benchmark in Table II. 



those cases where only an asymmetry is provided in [7], 
we multiply by the corresponding LO SM cross section 
to obtain A<r SM . For \i = 4m t , the A<j sm, s are given by 
the sum of their NLO values, obtained using MCFM [62], 
and estimates of their EW corrections, obtained by ex- 
trapolating the ones given in [7] . 

We exhibit the results of our scans in the (mz>,g u t) 
plane in Figs. 3-5, for three renormalization scales 



ji = m t /2^m t ^2m t . For each scale choice, results are 
shown for three representative branching ratios, Br(Z' — >• 
iu) = 1,1/3,1/4. For each scan (fixed /i and Br), the 
minimum x 2 and the corresponding values of (mz>,g u t) 
are indicated. In the plots, the la (2a) region of coverage 
probability satisfies X 2- Xmin < 2.3 (5.99), corresponding 
to variation of 2 parameters, mz> and g ut . The plots also 
show the regions which are compatible with the recent 
ATLAS bounds on top+jet resonance production [34], to 
be discussed in the next section. 
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Figure 5: Same as Fig. 3 for \i — 2m t , with the yellow crosses 
corresponding to the higher mz f benchmark in Table II. 



The best fit points tend to lie at low m^/, near 225 
GeV. Table II lists our best fit points, together with a few 
examples at heavier mz> near 300 GeV, which are marked 
with yellow crosses in the plots. The table also includes 
the best fit points for \i = Am t and Br = 1/4, 1/3 and a 
heavier example at this scale for Br = 1/4 (the Xmin = 3.1 
and 3.3 for Br = 1/4 and 1/3, respectively, and \ 2 — 5-4 
for the heavier example). For comparison, for the SM 
we find the x 2 values 12.1 (/i = rat/2), 12.9 (rat), 14.4 
(2rat), 15 (4rat). Thus, we find that the Z r model can 



significantly improve the fit. In particular, it can lower 
the x 2 by up to « 5, 8, 10, 12 for /i = rat/2, rat, 2rat, 4ra t . 

Regardless of scale choice, we find that the % 2 is mini- 
mized for Br(Z' — >• tu) « 1/4. The quality of the fits im- 
proves as the renormalization scale is increased, as does 
the overlap with regions that are consistent with the AT- 
LAS tj constraint. To first approximation, this is be- 
cause as the scale is increased: i) the LO NP quantities 
Aa tt,z and Aa z t decrease as one power of a s: while 
the LO SM cross sections in the denominators of the 
A' FB s decrease like a 2 s1 leading to an overall increase in 
the forward-backward asymmetries. A byproduct is that 
at larger /i the observed asymmetries can be reproduced 
with smaller couplings; ii) the magnitudes of the LO NP 

cross sections o\^ t Z and a z J: decrease as one power of 
a s which, when combined with the possibility of smaller 
couplings, makes it easier to satisfy the total cross sec- 
tion bounds. For example, for Br(Z' — >> tu) « 1/4 and 
/i = rat the minimum y 2 per d.o.f. (pi in Table II) cor- 
responds to a reasonable p- value of 17% or 1.4a, whereas 
for /i = 2rat it improves to 27% or 1.1 a (ps in Table II). 
Note that we have 3 d.o.f., as we fit 6 observables with 3 
parameters, i.e., the coupling, the mass and the branch- 
ing ratio. Reasonable fits are also possible for heavier Z' 
masses. For example, for mz> = 276 GeV and \i = 2ra t , 
Br = 1/3 (ps in Table II) the p-value is 23% or 1.2<j, and 
the ATLAS tj bound is satisfied. 

In Fig. 6 we compare the do j dm t i spectra for the 
SM and the five benchmarks listed in Table II 

with the measured Tevatron [11] and LHC [25] spectra. 
For the LHC the relative spectrum, (1/cr) da/dm t ^ is 
shown. The SM predictions are NLO for the Tevatron, 
and NLO+NNLL for the LHC [56]. The differences be- 
tween the SM predictions and the four benchmarks are 
smaller than the experimental errors (grey bands), with 
the exception of the last bin in the LHC m t i spectrum. 
The SM exceeds the ATLAS collaboration measurement 
by rsj la in this bin, and the five benchmark points ex- 
ceed it by ~ 1.5a (whether the "deficit" relative to the 
SM is a genuine effect, or is just the result of the last 
bin typically being statistically challenged will become 
clearer once the 8 TeV data is analyzed). Thus, the LHC 
is beginning to probe the high mass tail that results from 
exchange of the Z' in the t-channel [57]. 

We end this section with our most important result, an 
illustration in Fig. 7 of the crucial role played by associ- 
ated Z' production in reconciling the Tevatron and LHC 
ti asymmetries. The red (black) scatter points in the 
(AfbiAc) and (A^|> h , Ac) planes correspond to the la 
(2cr ) regions of coverage probability in the correspond- 
ing plots of Figs. 3-5. In addition we show in blue, 
whenever possible, the subset of points in the la regions 
that are also compatible with the ATLAS top+jet res- 
onance bound (see next section for details). The green 
points, obtained from the red and blue ones by setting 
Br(Z 7 — )> tu) = 0, clearly illustrate the sizable positive 
shifts in Ac which occur when associated Z' production 
is turned off. Moreover, while the red/blue la scatter 
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a(pp - 


* z't) [pb] 


a(pp — >> 


Z't) [pb] 




TEV 


LHC 


TEV 


LHC 


Pi 


0.57 


61.6 


0.57 


5.74 


P2 


0.19 


40.2 


0.19 


3.13 


P3 


0.39 


47.5 


0.39 


4.39 


P4 


0.079 


25.9 


0.079 


1.89 


P5 


0.24 


40.3 


0.24 


3.43 



Table III: The LO pp —> Z't and pp — »• Z't cross sections at 
the Tevatron (TEV) and the LHC for Table II benchmarks. 



making the analysis irrelevant for our case. 

The "wrong flavor" Z' — » ut decay could also give a 
contribution to the same-sign top pair cross section from 
pp —> Z't production. By our assumption Z' only decays 
into ui, as is the case in the flavor symmetric models, so 
this is not a problem. The experimental bounds from [48] 
are saturated for Br(Z' — » ut) < 3 — 5%, depending on 
the benchmark (the bounds in [49] and [50] are even less 
severe). They thus need to be an order of magnitude 
smaller than the allowed decay channel, Br(Z' — >• ut) ~ 
25%. 



B. Jet multiplicities 



points tend to lie well inside the experimental la band 
for Ac, the green points tend to lie above it. 



III. CONSTRAINTS FROM TOP+JET 
RESONANCE SEARCHES, AND JET 
MULTIPLICITY DISTRIBUTIONS 

A. top+jet resonance 

The Z' models we discuss can be searched for by hunt- 
ing for a bump in the f+jet distribution in the tij final 
state [35] . Such searches have been performed both at the 
LHC [33, 34] and the Tevatron [32]. The ATLAS search 
obtains a bound of cr(pp — >• tZ') + <j(pp — >> iZ') < 23, 14, 7 
pb for m z > = 200,300,400 GeV respectively. For a 
comparison we list in Table III the two cross sections, 
<j(pp — » tZ') and a(pp — >• iZ'), for the four benchmark 
points. We see that for Br(Z' — >> ut) = 1/4 the bench- 
mark points are not excluded. The comparison of the 
ATLAS bound with the regions in the g ut: mz> parame- 
ters space preferred by the tt asymmetry and cross sec- 
tion measurements are also shown in Figs. 3-5. The 
CDF search [32] is slightly less sensitive so that for the 
four benchmarks no bound on Br(Z' — »• ut) is obtained. 
The hard cuts in the CMS analysis [33] lead to a loss 
of sensitivity for light Z n s with masses below 400 GeV, 



Associated t + Z' — >• tij production would lead to some 
deviation in the jet multiplicity distribution for tt events 
relative to the SM prediction [37]. While it would be 
negligible at the Tevatron (cf. Table III), some sensitiv- 
ity to such deviations might be expected at the LHC. In 
Ref. [36] such a study was made for tt semileptonic candi- 
date events characterized by the number of 6-tags, lepton 
flavor and the number of jets. We focus on the double 
b- tagged subsample since it is subject to smaller back- 
ground contamination. The subsample is binned in terms 
of the number of jets. In the experimental analysis [36] 
a fit to the theoretical prediction was made, allowing the 
different backgrounds as well as the SM tt contributions 
to float in their normalization. Following this procedure, 
reasonable agreement with the SM Monte Carlo predic- 
tion was observed. The statistical error in bins of differ- 
ent jet multiplicities is of order a few percent. However, 
the background contamination is O(30%) or more in the 
bins with three or fewer jets and is O(10%) for the bins 
with four jets, and five or more jets. 

In order to check the viability of our mechanism we 
have performed the following semi-quantitative anal- 
ysis. Using MadGraph5 [51], Pythia6 .425 [52] and 
Fast Jet [53] we have produced an MLM matched ti+tij 
sample with and without the new physics contributions, 
mimicking the jet reconstruction and cuts employed by 
CMS. The computed jet multiplicities are compared with 
the experimental values obtained from Ref. [36] for the 
lepton plus jets channel (third column of Table 2). The 
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shown in blue. The green points are obtained from the la points by setting Br = 0. A yellow circle corresponds to the Xmin 
point, and its shift upon setting Br = 0. The yellow crosses correspond to a second heavier benchmark point, if indicated in 
Figs. 4, 5, and its shift for Br = 0. The gray bands are the experimental la intervals. 
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pp — »• tt + ttj 



SM 



0.12 



SM + Z'(pi) 
SM + Z'(p 2 ) 
SM + Z'(p 3 ) 
SM + Z x (p 4 ) 
SM + Z'(p 5 ) 



0.13 
0.14 
0.13 
0.14 
0.14 



0.22 
0.31 
0.25 
0.36 
0.31 



Table IV: Fraction of MLM matched ti+tij events with num- 
ber of jets equal to 5 or more for the SM and SM + Z' for the 
Table II benchmark points pi, ..,ps. 



results for the SM and for the five benchmark points in 
Table II are shown in Fig. 8. In none of the bins is the dif- 
ference between the SM only and the SM + NP samples 
larger than a few precent. . 

Of special interest for t + Z' — » ttj searches are the bins 
with large jet multiplicities. We thus analyze in more de- 
tail the fraction of events with 5 or more jets, with the 
results collected in Table IV. The fraction of events that 
have more than 5 jets is 0.12 for the SM, and is between 
0.13 and 0.14 for the five Table II benchmarks pi, ..,^5, 
well within the experimental errors. Isolating only the 
events with Z' associated production, the fraction is sig- 
nificantly larger, between 0.22 and 0.36 depending on 




Figure 9: Feynman diagram for the Z' contribution to 
the weak nuclear charge via the modification of the effective 
urZur coupling. 



the benchmark point. Because the pp — >• Z't — » tiX 
cross sections are substantially smaller than the total tt 
cross section, one sees that the impact on the combined 
SM+NP jet multiplicity is small. 



IV. CONTRAINTS FROM ATOMIC PARITY 
VIOLATION 

New massive bosons with parity violating couplings to 
light quarks can modify at one-loop the nuclear weak 
charge (Qw) measured in atomic parity violation (APV) 
experiments. Most constraining are the precise measure- 
ments with cesium ( 133 Cs) [63], which are in reasonable 
agreement with the most recent SM predictions [55] 



Q^(Cs) = 0.65(43). 



(14) 



This result can be interpreted in terms of constraints on 
the effective couplings of the SM Z boson to the light u 
and d quarks. It was argued in [54] that this leads to 
severe constraints on the Z r explanation of the Afb puz- 
zle. The argument is based on a calculation performed in 
a toy SU(2)x model [17], which contains two new spin-1 
states coupling to right-handed up-type quarks: Z' (V 
in the notation of [54]) with mostly flavor violating (u — t) 
couplings, and Z" (Z' in the notation of [54]) with mostly 
flavor diagonal (u — u and i — t) couplings. In this setup 
only the Z' is directly relevant for Afb- On the other 
hand, the strong APV constraints on the model are dom- 
inated by the Z" contributions. 

To demonstrate this we evaluate the modification of 
Qw in the minimal model (8). This can be viewed as the 
limit of the SU(2)x model with Z" heavy enough that 
it decouples (another possibility is to enlarge the sym- 
metry group to include a £7(1) factor, see below). We 
denote by Aclf>(u) the shift in the effective UrZur cou- 
pling induced by the Z' loop, cf. Fig. 9. Since the setup 
is non-renormalizable we impose a hard cut-off scale, A. 
This very roughly corresponds to the mass scale of the 
new degrees of freedom that will regulate the loop di- 
vergences found in the effective theory. In this way we 
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Figure 10: The Z' induced Aclr(u) contributing to APV in 
the minimal effective non-renormalizable model, as a function 
of the effective theory cut-off scale A for several representative 
Z' masses. 
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Figure 11: The Z' induced Aclr(u) contributing to APV 
in the renormalizable model, as a function of the vector-like 
quark mass m t > for several representative Z' masses (see text 
for details). 



obtain [54] 



Aa R (u) 



The dots denote additional Z" contributions, while 



-\9v 



mr 



16tt 2 



[FK 2 /m|0 + |log (A 2 /™?)] 



(15) 



where F{x) = (x— 1 — loga;)/(l — x) 2 . The resulting shift 
to the nuclear weak charge is given by 



AQ 



w 



-2(2Z n + N n )Aa R (u), 



(16) 



where Z p and N n are the number of protons and neutrons 
in a nucleus, respectively. One can thus translate Eq. (14) 
to a 1-cr preferred region for Aclr(u) = —0.17(11) x 10 -2 . 
This is to be compared with the predictions in the min- 
imal model as shown in Fig. 10. We have factored out 
the trivial dependence on the coupling g ut (note that in 
our best-fit regions \g u t\ ^ 1) an d plot Aa R (u)/\g ut \ 2 as 
a function of the effective theory cut-off A for several rep- 
resentative Z' masses. We see that for mz> > 300 GeV 
the APV bound can be satisfied even for cut-off scales 
above 1 TeV. 

The conclusions drawn in the effective non- 
renormalizable model can be further verified in a 
renormalizable ultraviolet completion of the SU(2)x 
model [54], where the Z' (and Z") loop divergences 
are cut-off by contributions of an additional heavy 
vector-like quark (t') with the SM gauge quantum 
numbers of the right-handed top. This gives 



Aa R (u) 



-\9u 



16tt 2 



-F 1 {m 2 t /m 2 zl 



m 2 ,/m 2 z ,) + ... 



(17) 



Fi(x,y) 



1 6 — 3x — 3y 

4 + (l-x)(l-y) 

(x 2 — 2x + 4) log x (2x 2 — 8x) log x 



(1-x) 2 
-2y + 4) logy 



(y 2 



(1- 

(2y 2 



x)(x - y) 
- 8y) logy 



(i-y) 2 



(l-y)(y-x) 



(18) 



In the m t > ^> mz> limit the above expression reproduces 
the log A dependence in Eq. (15) with A = m t ' . Again 
considering this as an illustrative setup (equivalent to 
the full SU(2)x model with a decoupling flavor diagonal 
Z") we plot Aclr(u) /\g ut \ 2 as a function of the t' mass 
(m t >) for several representative Z' masses in Fig. 11. We 
observe that the Z' contributions in the renormalizable 
model are even smaller than estimated in the effective 
theory and that the current APV constraints can be per- 
fectly satisfied even for mz> — 200 GeV (where the Afb 
fit prefers \g ut \ ~ 0.6) and t' masses well above 1 TeV. 

Finally, let us briefly discuss the viability of decou- 
pling Z" contributions to Aa^(w) in the SU(2)x model 
used above. The first approach to obtain this limit is to 
make the Z" heavy. However, the required hierarchical 
SU(2)x breaking pattern requires large scalar represen- 
tations [22]. Another possibility is to impose an OZI-like 
alignment of flavor diagonal vector mass eigenstates such 
that they couple predominantly to the u — uoii—t cur- 
rent and not both simultaneously, thus suppressing the 
flavor diagonal contributions to Aa^(w). For example, 
one can enlarge the symmetry group to include extra 
U{\) factors (i.e. to SU(2) X x U(l) x or U(2) x ). An 
explicit realization of the OZI alignment is then possible 
already with two scalar doublets breaking the symmetry 
group completely, by judiciously choosing their U(l) x 
charges. The desired OZI alignment would also follow 
dynamically from u uj — 0" mixing in a strong interaction 
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realization of the Z' model [38] in which, as previously 
noted, the Z' is the composite analog of a K* meson re- 
siding in a flavor nonet. The mixing of the [it] flavored 
"</>" and the [uu + cc] flavored "a;" would be suppressed, 
as in QCD. 

V. LEPTON BASED ASYMMETRIES 

An interesting Tevatron observable that is sensitive to 
the chiral nature of the NP model is the lepton based 
asymmetry 

A °i+(Vi+ > fy-<ri-(Vi- > 0) nq x 
l ~ a l+ (r ]l+ >0) + a l -(r ]l ->0y [ ] 

where a t ± (rji± > 0) are the tt cross sections for semilep- 
tonic t decays in which the / ± lepton pseudorapidity is 
positive. Similarly, the di- lepton asymmetries can be de- 
fined. The inclusive asymmetry has been measured to 
be A t = (15.2 ± 4.0)% by the DO collaboration [3], and 
A x = (6.6±2.5)% by the CDF collaboration [2]. The CDF 
collaboration has also measured the asymmetry in two 
m t t bins, obtaining Ai(m t t < 450 GeV) = (3.7 ± 3.1)% 
and Ai(m t t > 450 GeV) = (11.6±4.2)%. The SM predic- 
tions are Af M = (3.6 ± 0.2)%, Af M (m t t > 450 GeV) = 
(6.4 ± 0.5)% and Af M (m t i < 450 GeV) = (1.9 ± 0.1)% 
[7]- 

The lepton based asymmetry is an important observ- 
able because it gives an independent probe of the NP 
model explanations for the A^b measurements. In or- 
der to understand the physics that determines A\ in a 
qualitative manner it is best to consider two kinemati- 
cal regions for the tt invariant mass: near threshold, and 
well above threshold. Near threshold A\ is predominantly 
controlled by the chirality of the initial state quarks [67] . 
Well above the threshold it is determined by the polar- 
izations of the tops [68-70] . In our framework the Z' cou- 
ples to right-handed quarks. Thus, we expect a negative 
value for A\ near threshold, from the interference between 
the SM amplitude and the NP ^-channel amplitude. The 
contribution to A\ from associated Z' production is small 
due to its small cross section at the Tevatron. At high 
m t t we expect a positive value of A\ because more right- 
handed tops are produced from the NP interactions. This 
implies that there will be some cancelation between these 
two kinematical regions in the inclusive A\. Therefore, 
the exact value of the inclusive asymmetry will be rela- 
tively sensitive to the details of the NP parameter space. 

Numerically, employing the prescription in (10), we 
find that the inclusive values of A\ lie between 4% to 6% 
for all of our benchmark points (pi, ..,^5), consistent with 
the CDF measurement and somewhat lower than the DO 
measurement. We also find that A\(m t i < 450 GeV) ~ 
-7 to -8% and Ai(m t t > 450 GeV) - 20 - 23% for all 



of our benchmark points (pi, ..,^5). Thus, the tendency 
of A\ to increase with m t i is indeed present in the Z' 
models, but the rise is more pronounced than observed 
by CDF. 

VI. DISCUSSION 

In conclusion, we have shown that associated produc- 
tion of a vector state, pp — » tV, where V — » ij, gives a 
negative contribution to the tt charge asymmetry Ac at 
the LHC, while it has a negligible effect on the forward- 
backward asymmetry Afb measured at the Tevatron. 
As is well known, the same underlying interaction also 
gives rise to positive contributions to Afb and Ac, via 
exchange of the vector in the t-channel. Thus, Ac can 
have SM like values, while Afb can be anomalously large 
and positive, in agreement with experiment. We have dis- 
cussed this quantitatively for the example of a Z', where 
the Z' decays to ut roughly 25% of the time. The main 
result is summarized in Fig. 7, showing a substantial neg- 
ative shift in Ac and a much smaller shift in Ap b •> bring- 
ing the Z 1 models into agreement with their measured 
values. 

Finally, we have shown that the Z* model is not ruled 
out by other searches: the top+jet resonance search in 
the tij final state, the jet multiplicity distributions in tt 
production, and atomic parity violation measurements. 
It would be interesting to search for the dominant Z' 
decay channel candidates, e.g., invisible final states, or 
t-\-n jet final states. A test of the t— channel Z' exchanges 
is provided by lepton based asymmetries. 

Note added: While our work was being submitted to 
the arXiv, similar conclusions appeared in Ref. [71]. 
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